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ARTICLE INFO ABSTRACT 
Article history: This paper reviews the development of latent heat thermal energy storage systems studied detailing 
Received 4 August 2009 various phase change materials (PCMs) investigated over the last three decades, the heat transfer and 


Accepted 19 October 2009 enhancement techniques employed in PCMs to effectively charge and discharge latent heat energy and 


the formulation of the phase change problem. It also examines the geometry and configurations of PCM 
Keywords: containers and a series of numerical and experimental tests undertaken to assess the effects of 
Phase change material (P CM) | parameters such as the inlet temperature and the mass flow rate of the heat transfer fluid (HTF). It is 
isn peli Amie energy storage (LHTES) concluded that most of the phase change problems have been carried out at temperature ranges between 
Conduction/convection heat transfer 0 °C and 60 °C suitable for domestic heating applications. In terms of problem formulation, the common 
approach has been the use of enthalpy formulation. Heat transfer in the phase change problem was 

previously formulated using pure conduction approach but the problem has moved to a different level of 
complexity with added convection in the melt being accounted for. There is no standard method (such as 
British Standards or EU standards) developed to test for PCMs, making it difficult for comparison to be 
made to assess the suitability of PCMs to particular applications. A unified platform such as British 
Standards, EU standards needs to be developed to ensure same or similar procedure and analysis 
(performance curves) to allow comparison and knowledge gained from one test to be applied to another. 
© 2009 Elsevier Ltd. All rights reserved. 
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1. Introduction 
The quest for new technologies to avert the growing concern 
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E-mail address: fb.agyenim@ulster.ac.uk (F. Agyenim). the high cost of energy and new power plants has been a scientific 
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Nomenclature E emissivity, emmittance 
| o . T dimensionless time 
A area (m°); thermal diffusivity ratio, A = as/ar 
a. equation coefficient Subscripts and superscripts 
Bi Biot number, Bi = hr,/À C convective 
B solid subcooling parameter, B = A;AT;/A FAT, Cd conductive 
b width of a rectangular PCM container (m) P cylindrical 
C instantaneous height-averaged cooling rate ie bel 
G specific heat capacity (kJ/kg K) e ne 
d diameter (m) eff effective 
Fo Fourier number eq equivalent 
f liquid fraction of PCM melt f fluid 
g acceleration due to gravity (m/s) fiber fiber 
Gr Grashof number fin fin 
al enthalpy (J) fric friction 
h heat transfer coefficient (W/(m? K)) H heat transfer fluid 
He instantaneous, height-averaged melting rate i inner (inlet) 
HTF heat transfer fluid l liquid 
k thermal conductivity (W/m) lag lag 
l length (m) i melting 
M instantaneous, height-averaged melting rate max maximun 
m mass (kg) o outer (outlet) 
mmass flow rate (kg/s) PCM 
kf Prandtl number : r rectangular container 
PCM phase change material ref reference 
Q energy or heat flow (J) ‘ solid 
q heat flux (W/m?) i tube 
R thermal resistance (K/W) 
Ra Rayleigh number 
= nnn a concern over the last three decades. Central to the problem is the 
. . need to store excess energy that would otherwise be wasted and 
Ste Stefan number or liquid superheating, Ste = CAT;/ also to bridge the gap between energy generation and consump- 
AT, tion. Latent heat thermal energy storage is particularly attractive 
T temperature (°C or K) technique because it provides a high energy storage density. When 
AT temperature difference (°C or K) compared to a conventional sensible heat energy storage systems, 
t time (s) latent heat energy storage system requires a smaller weight and 
At time step (s) volume of material for a given amount of energy. In addition latent 
v velocity components (m/s) heat storage has the capacity to store heat of fusion at a constant or 
V volume (m?) near constant temperature which correspond to the phase 
7 velocity (m/s) transition temperature of the phase change material (PCM). The 
x lisci volume trictionof fibers study of phase change materials was pioneered by Telkes and 
g i . : Raymond [1] in the 1940s but did not receive much attention until 
XYZ KY z-component in a Cartesian coordinate system the energy crisis of late 1970s and early 1980s where it was 
Ax space increment (m) extensively researched for use in different applications especially 
wW dimensionless geometric parameter defined by 1 + for solar heating systems [2-33]. Less attention was however given 
w = A/V to phase change storage research after the energy crisis until it 
Z fin height (m) emerged again recently. Since the late 70s a number of studies 


Greek symbols 


ô 


FAT ErD wR 


fin thickness (m); dimensionless front radius 
(interface position) 

thermal diffusivity of PCM (m2/s) 
coefficient of volumetric expansion (K~!) 
density (kg/m?) 

latent heat of fusion (kJ/kg) 

angle (°) 

dimensionless front radius, ô = rfr, 
dimensionless temperature difference 
dynamic viscosity (m?/s) 

viscosity (kg/(m s)) 


have been conducted to assess the overall thermal behaviours of 
latent heat thermal storage systems [34-59]. Studies of phase 
change systems have investigated design fundamentals, system 
and process optimization, transient behaviour, and field perfor- 
mance. The research and development has been broad based and 
productive, concentrating on both the resolution of specific phase 
change materials and problems and the study of the characteristics 
of new materials. The major disadvantage, as reported by many 
researchers has been the low thermal conductivities possessed by 
many PCMs, leading to low charging and discharging rates 
(especially for the organic based materials). The development of 
a latent heat thermal energy storage system therefore involves the 
understanding of heat transfers/exchanges in the PCMs when they 
undergo solid-to-liquid phase transition in the required operating 
temperature range, the design of the container for holding the PCM 
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and formulation of the phase change problem. This paper reviews 
the studies conducted on phase change materials for different 
applications and reported in the open literature. 


2. Phase change materials (PCMs) 


Energy storage may be in the form of sensible heat in a liquid or 
solid medium, as heat of fusion (latent heat), or as chemical energy 
or products in a reversible chemical reaction. The classification of 
energy storage and the materials used are detailed in [51]. 
Chemical energy storage has not as yet been used in practical 
applications and both technical and economical questions have yet 
to be answered for some of the possibilities proposed. To date most 
of the studies conducted on storage materials have concentrated 
on sensible and latent heat storage systems. Studies conducted to 
compare phase change and sensible heat storages have shown that 
a significant reduction in storage volume can be achieved using 
PCM compared to sensible heat storage. Results of studies from 
Morrison and Abdel-Khalik [48] and Ghoneim [60] show that to 
store the same amount of energy from a unit collector area, rock 
(sensible heat storage material) requires more than seven times 
the storage mass of Paraffin 116 Wax (P116-Wax), five times the 
storage mass of medicinal paraffin and more than eight times the 
storage mass of Na2SO,-10H20. 


2.1. PCMs investigated 


Several authors have carried out investigation into a wide range 
of PCMs, subdividing them into organic, inorganic, eutectics PCMs 
[61]. The main criteria that govern the selection of phase change 
heat storage materials are [51]: 


e Possess a melting point in the desired operating temperature 
range (temperature range of application). 

e Possess high latent heat of fusion per unit mass, so that a smaller 
amount of material stores a given amount of energy. 

e High specific heat to provide additional significant sensible heat 
storage effects. 

e High thermal conductivity, so that the temperature gradients for 
charging and discharging the storage material are small. 

e Small volume changes during phase transition, so that a simple 

container and heat exchanger geometry can be used. 

Exhibit little or no subcooling during freezing. 

Possess chemical stability, no chemical decomposition and 

corrosion resistance to construction materials. 

e Contain non-poisonous, non-flammable and non-explosive 
elements/compounds. 

e Available in large quantities at low cost. 


Table 1 lists the current companies that commercially produce over 
100 PCMs. In addition to these, several PCMs have been proposed or 
studied by different researchers. A detailed list of PCMs studied or 
proposed for study can be found in [51,62-65]. Table 2 gives the 


Table 1 

Commercial PCM manufacturers in the world. 
Manufacturer PCM temperature Number of 

range PCMs listed 

RUBITHERM (www.rubitherm.de) —3 °C to 100°C 29 
Cristopia (www.cristopia.com) —33 °C to 27°C 12 
TEAP (www.teappcm.com) —50°C to 78°C 22 
Doerken (www.doerken.de) —22°C to 28°C 2) 
Mitsubishi Chemical (www.mfc.co.jp) 9.5°C to 118°C 6 
Climator (www.climator.com) —18°C to 70°C 9 


EPS Ltd (epsltd.co.uk) —114°C to 164°C 61 


Merck 


thermophysical data for a range of typical and important PCMs 
studied or proposed for study by previous authors. 


2.2. Classification of PCMs to melting temperature range and 
application area 


The selection of an appropriate PCM for any application 
requires the PCM to have melting temperature within the practical 
range of application. Several application areas have been proposed 
for PCMs studied. Table 2 lists some of the target application areas 
for the selection of PCMs for study, as provided by the references in 
the open literature. From Table 2, it can be seen that most of the 
research on phase change problems have been carried out within 
the temperature range 0-65 °C suitable for domestic heating/ 
cooling (Table 3). 


3. PCM containers 


Once the PCM has been selected based primarily on the 
temperature range of application, the next most important factors 
to consider are: 


(i) the geometry of the PCM container and 

(ii) the thermal and geometric parameters of the container 
required for a given amount of PCM. 

Each of these factors has a direct influence on the heat transfer 

characteristics in the PCM and ultimately affects the melt time and 

the performance of the PCM storage unit. 

To ensure long-term thermal performance of any PCM system, 
the size and shape of the PCM container must correspond to the 
melting time of the PCM and the daily insolation at a given 
location, if the source of energy is a solar collector. PCMs are 
typically placed in long thin heat pipes [102], cylindrical containers 
[100,103] or rectangular containers [66,91]. A survey of previously 
published papers dealing with LHTES reveals that two geometries 
commonly employed as PCM containers are the rectangular and 
cylindrical containers. The most intensely analysed LHTES unit is 
the shell and tube system, accounting for more than 70%. This is 
probably due to the fact that most engineering systems employ 
cylindrical pipes and also heat loss from the shell and tube system 
is minimal. Fig. 1 gives the schematics of the cylindrical and 
rectangular containers. Table 4 lists some of the references to the 
studies undertaken employing the two most common container 
geometries. 


3.1. Configurations of cylindrical PCM containers 


Three modes of cylindrical PCM container configurations are 
distinguished. The first is where the PCM fills the shell and the heat 
transfer fluid flows through a single tube (Fig. 1a) [60,80,98-100] 
designated the pipe model. In the second model the PCM fills the 
tube and the HTF flows parallel to the tube (Fig. 1b) [56,107,80]. 
According to Esen et al. [80], who studied the two models 
theoretically by comparing the effects of various thermal and 
geometric parameters; cylinder radii, total PCM volume, mass flow 
rates and inlet temperatures of HTF on the storage time 
recommended the pipe model because it recorded a shorter melt 
time. This was because the thicker the PCM mass, the longer the 
melt time of the PCM. An added advantage for the pipe model not 
mentioned is the fact that the pipe model has a lower heat loss rate 
to the environment because most heat supplied from the center 
ends up heating the PCM. The third cylinder model is the shell and 
tube system [60,99] commonly used to improve heat transfer in 
PCMs. Agyenim et al. [99] conducted an experimental energy 
storage system to compare horizontal shell and tube heat 
exchanger (4 tubes) and a pipe model incorporating a medium 


Table 2 


Thermophysical properties of PCMs investigated for different applications. 


Compound Reference(s) Melting temp, Tm (°C) Heat of fusion, Specific heat capacity, Thermal conductivity, Density, p (kgm) 

À (kj kg7!) G (kJ kg 'K !) k (Wm™! K!) 
Water-ice [66-70] 0 335 42 2.4 (liquid) 0.6 1000 
GR25 [74] 23.2-24.1 45.3 1.2 (solid) 1.2 (liquid) - - 
RT25-RT30 [72,73] 26.6 232.0 1.80 (liquid) 1.41 (solid) 0.18 (liquid) 0.19 (solid) 749 (liquid) 785 (solid) 
n-Octadecane [55,56,74-79] 27.7 243.5 2.66 (liquid) 2.14 0.148 (liquid) 0.190 (solid) 785 (liquid) 865 (solid) 
CaCl2-6H20 [66,80,81] 29.9 187 2.2 (liquid) 1.4 (solid) 0.53 (liquid) 1.09 (solid) 1530 (liquid) 1710 (solid 
Rubitherm RT [37] 
Na2SO,-10H20 [60,80,81] 32, 39 180 2.0 (liquid) 2.0 (solid) 0.15 (liquid) 0.3(solid) 1460 (solid) 
Paraffin wax [72,80,81] 32-32.1 251 1.92(solid) 3.26(liquid) 0.514 (solid) 0.224 (liquid) 830 
Capric acid [82,83] 32 152.7 - 0.153 (liquid) 878 (liquid) 1004 (solid 
Polyethelene glycol 900 (PEG900) [84] 34 150.5 2.26 (liquid) 2.26 (solid) 0.188 (liquid) 0.188 (solid) 1100 (liquid) 1200 (solid 
Lauric-palmitic acid (69:31) eutectic [85] 35.2 166.3 2.41 (liquid) 1.77 (solid) 
Lauric acid [82,85-87] 41-43 211.6 2.27(liquid) 1.76(solid) 1.6 1.76(solid) 0.862 (liquid) 
Stearic acid [61,87-89] 41-43 (67-69)* 211.6 2.27 (liquid) 1.76 (solid) 1.60 (solid) 862 (liquid) 1007 (solid 
Medicinal paraffin [60] 40-44 146 2.3 (liquid) 2.2 (solid) 2.1 (liquid) 0.5 (solid) 830 (solid) 
Paraffin wax [90] 40-53 
P116-Wax [55,60,80,81] 46.7-50 209 2.89 (liquid) 2.89 (solid) 0.277 (liquid) 0.140 (solid) 786 (solid) 
Merck P56-58 [53] 48.86-58.06 250 2.37 (liquid) 1.84 (solid) - - 
Commercial paraffin wax [94] 52.1 243.5 - 0.15 809.5 (solid) 771 (liquid) 
Myristic acid [35,58,95] 522 182.6 - - - 
Paraffin RT60/RT58 [93,94 55 to 60 214.4-232 0.9 0.2 775 (liquid) 850 (solid) 
Palmitic acid [83,85,87,89,92,95] 57.8-61.8 185.4 0.162 (liquid) 850 (liquid) 989 (solid) 
Mg(NO3)2-6H20 [71,88 89 162.8 - 0.490 (liquid) 0.611 (solid) 1550 (liquid) 1636 (solid) 
RT100 [96] 99 168 2.4 (liquid) 1.8 (solid) 0.2 (liquid) 0.2 (solid) 770 (liquid) 940 (solid) 
MgCl,-6H20 [71,97 116.7 168.6 2.61 (liquid) 2.25 (solid) 0.570 (liquid) 0.704 (solid) 1450 (liquid) 1570 (solid) 
Erythritol [50,98-100] 1177 339.8 2.61 (liquid) 2.25(solid) 0.326 (liquid) 0.733 (solid) 1300 (liquid) 1480 (solid) 
Na/K/NO3 (0.5/0.5) [54,65 220 100.7 1.35 0.56 1920 
ZnCl/KCI (0.319/0.681) [54,65 235 198 - 0.8 2480 
NaNO, [54,65 310 172 1.82 05 2260 
KNO3 [54,65 330 266 1.22 0.5 2110 
NaOH [54,65 318 165 2.08 0.92 2100 
KOH [54,65 380 149.7 1.47 0.5 2044 
ZnCl [54,65 280 75 0.74 0.5 2907 
LiF-CaF; (80.5:19.5) mixture [101] 767 816 1770 (liquid) 1.770 (liquid) 1.70 (liquid) 3.8 (solid) 2390 (liquid) 2390 (solid) 


- Not available. 


819 
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Table 3 
Target application areas for some PCMs studied. 


Temperature PCMs studied/melting temperature (°C) Target application area (rationale behind selection of PCM) 

range (°C) 

0-65 °C Paraffins (—3 to 64), water/ice/0, stearic acid/41-43, Storage for domestic heating/cooling (see list in Table 1). Passive 

n-octadecane/27.7 storage in bio-climatic building/architecture. Thermal storage of solar 

energy. Application in off-peak electricity for cooling and heating. 
Protection of electrical devices. 

80-120°C Erythritol/117.7; RT100 (99); MgCl:-6H20 (116.7) Storage for the hot-side of LiBr/H20 absorption cooling system with 
generator temperature requirements of less that 120°C [100]. 

>150 NaNO;/310, KNO3/330 NaOH/318, KOH/380 ZnCl,/280 Storage for solar power plants based on parabolic trough collectors 


and direct steam generation. 


(a) Pipe model 


(b) Cylinder model 


(c) Shell and tube model 


(d) Rectangular or slab container 


Fig. 1. Classification of commonly used PCM containers in terms of the geometry and configuration. 


temperature phase change material (erythritol) with a melting 
point of 117.7 °C. The thermal characteristics in the systems were 
analysed using isothermal contour plots and temperature time 
curves. Temperature gradients along the three directions of the 
two systems; axial, radial and angular directions were also 
analysed and compared. Heat transfer in the shell and tube 
system was found to be dominated by the effect of multiple 
convective heat transfer compared to conductive heat transfer in 
the pipe model. The temperature gradients recorded in the axial 
direction for both the pipe and shell and tube systems during the 
change of phase were reported to be 2.5% and 3.5% that of the radial 
direction respectively, indicating essentially a two-dimensional 
heat transfer in both systems. The onset of natural convection 
through the formation of multiple convective cells in the shell and 
tube system significantly altered the shape of the solid liquid 
interface fluid flow and indicated complete melt time within 5h 
compared to more than 8h for the pipe model. The authors 
recommended the shell and tube system for the charging of PCMs. 


Table 4 
References to container geometry employed in the study of PCMs. 


Geometry Reference(s) 


[60,66,74,76,78,80,81,86,93,95, 
97,102,106,107] 

Cylindrical (shell and tube) [56,75,90,101,104,105] 
Rectangular/slab [66,68,72,74,77,79,91,108] 
Spherical [74] 


Cylindrical (concentric annulus) 


3.2. Counter-current and parallel HTF flow directions 


In a cylindrical container assembly, two possibilities exist for 
the flow direction of the heat transfer fluid during charging and 
discharging of the PCM energy. The two modes are either the hot 
and cold fluids are introduced into the tube from the same end 
(parallel flow) or the hot and cold fluids are introduced from the 
opposite ends (counter-current flow) of the heat transfer tube 
during charging and discharging respectively. Fig. 2 illustrates the 
schematic diagram of the parallel and counter-current flow 
principles. For each pair, the upper arrow represents the direction 
of HTF flow during charging and the lower arrow represents 
discharge direction of the HTF. 

Results from the numerical simulations conducted by Gong and 
Mujumdar [101] to investigate the effect of the parallel and 
counter-current flow modes using a mixture of 80.5% LiF and 19.5% 
CaF2 as a PCM and He/Xe mixture as a HTF showed that parallel 
flow increases the energy charge/discharge rate by 5% more than 
counter-current flow. This was because the penetration depth of 
the solid-liquid phase change interface during the charge/ 
discharge was larger due to higher temperature difference at 
the fluid inlet if hot and cold fluids enter from the same end. In 
addition, supercooling of the PCM did not occur in the fluid inlet 
region and heat transfer between the heat transfer fluid and the 
PCM did not deteriorate. On the other hand, counter-current flow 
for the charge/discharge processes produced significant super- 
cooling of the PCM in the inlet region of the cold fluid. 
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Charging 


Discharging 


(a) Parallel HTF flow 


l 


Charging Charging 


Discharging Discharging 


(b)Counter-current HTF flow 


Fig. 2. The physical model illustrating parallel and counter-current HTF flow in a 
shell and tube system. 


3.3. Investigating parameters in PCM containers 


Apart from the container geometry and configuration, various 
thermal and geometric parameters are known to affect the thermal 
performance of LHTES. Table 5 lists the references to some of the 
operating parameters investigated and their key findings. 

Key findings from numerical predictions, tested against 
experimental data for various operating parameters have shown 
that, in order to optimize the performance of a PCM storage unit, 
the thermal and geometric parameters listed in Table 5 must be 
carefully selected. Literature reviewed has shown a superior 
performance using the shell and tube configuration followed by 
the pipe model with the PCM at the shell side and the heat transfer 
fluid flowing through the center. It has therefore been recom- 
mended by many authors including [60,75,99,118]. 


4. Heat transfer in PCMs and enhancement techniques 


Most PCMs have unacceptably low thermal conductivity, 
leading to slow charging and discharging rates, hence heat transfer 
enhancement techniques are required for most LHTES applica- 
tions. Several studies have been conducted to study heat transfer 
enhancement techniques in phase change materials (PCMs) and 
include finned tubes of different configurations [59,69,94,97,98, 
102,109-112], bubble agitation [112], insertion of a metal matrix 
into the PCM [113,114], using PCM dispersed with high 
conductivity particles [115], micro-encapsulation of the PCM 
[116,117] or shell and tube (multitubes) [99,118]. Fig. 3 shows the 
pictures and drawings of some of the common heat transfer 
enhancement techniques studied. Table 6 lists some examples 
where enhancement techniques have been investigated. 

Majority of the heat enhancement techniques have been based 
on the application of fins embedded in the phase change material. 
This is probably due to the simplicity, ease in fabrication and low 
cost of construction. This is followed by the impregnation of metal 


matrix into the PCM using high conductivity materials such as 
carbon fiber and brushes and multitubes. A selection of four studies 
where fins (circular and longitudinal), metal matrix and shell and 
tube systems approach have been adapted are described in Table 
7. It can be seen from Table 7 that different approaches were used 
in the analysis with different results and interpretations. 

The general observations drawn from the various studies 
demonstrate that, irrespective of the PCM used, the heat transfer 
characteristics of the PCMs can be improved using all of different 
enhancement techniques. Choi and Kim [97], Horbaniuc et al. 
[102], Velraj et al. [93] and Hamada et al. [75], employed 
different experimental setups, different container materials, and 
different PCMs to investigate the heat transfer enhancement 
characteristics of PCMs. 

In terms of performances of heat transfer enhancement 
techniques and systems used by Choi and Kim [97], Horbaniuc 
et al. [102], Velraj et al. [93] and Hamada et al. [75], the best 
enhancement technique as reported in the literature was that due 
to Velraj et al. [93] where the effective thermal conductivity 
calculated employing paraffin with lessing rings was ten times 
(2Wm-!K-1) greater than the thermal conductivity of paraffin 
(0.2 W m~! K~!). This must be qualified by the fact that different 
researchers employed different experimental setups, different 
container materials, and different PCMs to investigate the heat 
transfer enhancement characteristics (Table 8). For example 
during the charging and discharging stages of the experiment, the 
heat transfer fluid in the experiment from Horbaniuc et al. [102] 
passed over just a fraction of the whole setup rendering heat 
addition and removal inefficient. Results derived from any such 
setup cannot be considered to be optimal. 

Different researchers used different parameters to assess the 
heat transfer enhancement in the PCMs. Velraj et al. [93] evaluated 
the enhancement of the heat transfer using the effective thermal 
conductivity taken from a two-dimensional enthalpy-temperature 
governing equation which assumed no variation of temperature 
and thermal conductivity in the axial direction. Results were 
presented graphically using temperature-time curve and as such 
limits the application tending to other applications. Horbaniuc 
et al. [102] measured performances of fins in terms of the interface 
freezing stage and the time taken for complete solidification to be 
achieved using parabolic and exponential approximations. 
Hamada et al. [75] used the effective thermal conductivity 
proposed by Fukai et al. [126] to assess and compare results 
with the control system with no heat transfer enhancement. In the 
case of Choi and Kim [97], the key parameters used to assess the 
heat transfer enhancement of the circular finned system were the 
ratio of overall heat transfer coefficient in the finned and the 
unfinned tube systems. They reported a ratio of 3.5 for a surface 
area ratio of 3.2 between the finned and the unfinned tube systems. 
Acomparison was also made by deriving a relationship for the ratio 
of the total amounts of heat recovered in the finned and unfinned 
tube systems correlated with dimensionless parameters of Fourier, 
Stefan and Reynolds with a correlation coefficient of 0.994 and a 


Table 5 
Examples of common operating parameters investigated in the study of phase change containers. 
Reference(s) Parameter(s) investigated Findings 
[56] 1. PCM cylinder (shell) radii 1. Increasing mass flow rate increases the heat transfer rate and the molten volume fraction. 
2. Mass flow rate 2. Increasing inlet temperature rate increases the heat transfer rate and the molten volume fraction. 
3. HTF inlet temperature 
[66] Comparison of rectangular and For larger quantities of PCM (>2 kg), the melting time of cylindrical container was nearly twice 
cylindrical containers that of the rectangular one for equal volume and equal heat transfer. 
[74] Melt time Average and total heat fluxes (4852 kJ/m? and —4403 kJ/m?) approach the latent heat storage 
capacities per unit surface area as change in temperature approaches zero. 
[101] Effect of the parallel and Energy charge/discharge rate for parallel flow is 5% more than counter-current flow mode 


counter-current flow modes 
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(iii) Multitubes or (iv) Bubble agitation 
shell and tube 


(i) Longitudinal or 
axial fins 


(vi) Multitubes and (vii) Encapsulation (ix) Metal Matrix 
carbon brushes 


ý 


À : 19, í 
(x) Finned Re ctangular (xi) Graphite flakes (xii) Steel metal (xiv) Polyolefine 
Container ball capsules spherical balls 


(xv) Polypropylene (xvi) Module beam (xvii) PCM- (xviii) Compact flat 
flat panel Graphite panel 


Fig. 3. Heat transfer enhancement methods employed in phase change material research. 


Table 6 
Heat transfer enhancement techniques employed in PCMs. 
Reference(s) Nature of study Description of enhancement technique 
[119] Experimental LiH PCM encapsulated in smaller container and placed in a bigger cylindrical container 5% Ni added uniformly to LiH. 
[120] Analytical Aluminium and steel internal fins in the solidification of n-octadecane and salt-hydrate climsel. 
[118] Numerical and experimental Staggered cylindrical tubes (multitubes) in a shell and tube system using Indonesian traditional substance 
called Mikro as a PCM. 
[121] Numerical Plate type fins mounted on the vertical cooling surface of a rectangular vessel to study solidification of nitric 
molten salt. 
[112] Numerical and experimental Internal longitudinal fins to improve heat transfer in paraffin PCM in a cylindrical container. 
[122] Experimental A spiral twisted copper fins in a cylindrical tube to enhance heat transfer during the melting of stearic acid. 
[110] Numerical Melting of ice in a porous medium surrounded by finned surfaces. 
[79] Numerical and experimental System of 7 aluminium rectangular containers filled with n-octadecane and used as fins to improve 
heat transfer of PCM. 
[123] Analytical Simplified analytical model to study heat transfer of fins and honeycombs in PCMs. 
[124] Experimental Ultrasonic energy used to enhance heat transfer in Erythritol PCM. 
[125] Numerical Micro-encapsulation of phase change slurries to improve convective heat transfer in the PCM. 


[99,100] Experimental Internal longitudinal and circular fins and shell and tube system to improve heat transfer in Erythritol. 
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Table 7 
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A comparison of results of study using fins, metal matrix and shell and tube system. 


Reference(s) Heat enhancement Method Analysis and results Observation or recommendation 
technique 
[97] Circular fins Freezing of 0.36 kg weight Unfinned tube system:, 2 = 2.39x Radial temperature gradient in 
MgCl2-6H20 10-2 Fo°-765 Ste9528Re°195 Re < 7000 the finned tube system is uniformly 
larger than the unfinned tube system. 
Finned tube system: T =E The axial temperature difference for 
10-3 Fo°-781 Ste0-616 Re? 303" Re > 7000 the finned tube was smaller than for 
the unfinned tube in the liquid phase. 
[102] Longitudinal fins Solidification of PCM Parabolic approximation of temperature Ten fins would be sufficient to allow 
(Tm= 60°C) profile: t{a(t)]? = 225 a complete solidification of PCM with 
mf 2 Tm=60°C in an imposed freezin 
AT: -AT2 m p S 
2 Ieee ere E time of 180 min. 
2 fin 2Jalt) VAT fin 22/000) 
Exponential approximation of temperature 
profile: -Ø"4In® = a,xt-1/4, D = a and 
— __2Apsks 
ox 17,0 fin AT fin 
[93] Fins, rings and Solidification of Paraffin Recommended internal longitudinal fins and Effective thermal conductivity keg 
bubble agitation RT 60 and RT 58 lessing rings for solidification enhancement calculated by employing paraffin 
and the use of an evacuated tube for melting with lessing rings was 2Wm !°C-! 
when using paraffin. (ten times greater than the thermal 
conductivity of paraffin). 
[75] Carbon fiber Melting of PCM Estimation of effective thermal conductivities Heat brush/PCM composite was found 
and brushes n-octadecane using the semi-empirical equation from Fukai to be inferior to the fiber-chip/PCM 
k ite i f heat transfer 
et al.: “ell — 04561 -(1-X 2/3 composite in terms o 
4 { 4 fiver) } rate of the bulk but the former was 
(K fiber/Kp) + (1 — X fiber) found superior to the latter in terms 
of overall heat transfer rate. 
Table 8 


Container materials used by previous researchers for PCM heat transfer enhancement. 


Reference(s) PCM container material and position Heat transfer tube Enhancement technique 
[97] Pyrex glass/vertical Stainless steel tube 5 stainless steel circular fins 
[102] Horizontal - Longitudinal fins 

[93] Aluminium/vertical Aluminium tube Fins, rings, bubble agitation 
[75] Acrylic resin/vertical 4 steel tubes Carbon fiber and brushes 


standard deviation of 0.023 and 0.028, respectively. The above 
discussion illustrate the fact that there is no unified international 
or national standard methods (such as British Standards or EU 
standards) developed to test PCMs, making it difficult for 
comparison to be made to assess the suitability of PCMs to 
particular applications. A standard platform such as British 
Standards, EU standards needs to be developed to ensure same 
or similar procedure and analysis (performance curves) to allow 
comparison and knowledge gained from one test to be applied to 
another. Also contradictions exist in the thermophysical properties 
of PCMs provided especially the latent heat of fusion, thermal 
conductivity and densities in solid and liquid states. This is again 
due to the absence of unified certification standards and 
procedures. 


5. Characterising the effect of melting and solidification rates 
in PCMs 


To design PCM systems requires a good understanding of the 
fundamental heat transfer processes involved in accurately 
predicting the thermal performance of the PCM system. Based 
on both experimental and numerical investigations conducted to 
examine the thermal characteristics of LHTES systems, various 
correlations have been developed relating thermal performance 
and dimensionless numbers in given parametric domains. 
Dimensionless numbers are widely applied in modelling the 
PCM systems in order that the knowledge gained from one study 
can be extended beyond its source of acquisition. Among the 
thermal performance parameters for which correlations have been 


developed are melted volume fraction, temperature profile, melt 
time and melting rate. Table 9 gives the mathematical definitions 
of the common dimensionless numbers used in the study of LHTES 
systems and their significances to the phase change process. Table 
10 gives references to some of the correlations derived using 
dimensionless numbers. 

In spite of the significant literature available in quantifying 
melting and solidification rates, there has been no attempt to 
correlate available data. The two main reasons for this are: 


e individual authors used different phase change materials with 
different heat transfer characteristics. In the case where the same 
PCM has been used, the researchers employed different 
dimensionless parameter ranges and presentations making it 
difficult to cross-correlate between the characteristics influen- 
cing the heat transfer in specific PCMs. 

e some researchers presented results without using dimensionless 
parameters making it impossible to extend the knowledge to 
applications beyond the original source. 


For example the empirical correlation for the melted volume 
fraction from Wang et al. [84], for example predicted the 
experimental data with an average error of 9.3% and the deviation 
was attributed to the sensible heating that takes place at the very 
early times due to the presence of subcooling within the PCM. In 
the case of the time averaged Nusselt number, the correlation was 
found to fit the individual experimental data with an average error 
of 8.8%. The dependency of the melted volume fraction on Fo and 
Ste but minimal for Ra was that, the dimensionless groups Fo and 


F. Agyenim et al. / Renewable and Sustainable Energy Reviews 14 (2010) 615-628 623 


Significance (determination) 


Ratio of conductive to convective heat transfers resistance. 
Determines uniformity of temperature in solid. 

Ratio of the conductive thermal resistance to convective thermal 
resistance. Determines the ratio of actual heat transferred 

by a moving fluid to the heat transfer that would occur by conduction. 
Ratio of thermal capacity of the melted solid to the latent heat. 
Characterises heat flux into a body or system. 


Determines the onset of convection. Below a critical value, 
heat transfer is primarily conduction. 

Approximates the ratio of momentum diffusivity to thermal 
diffusivity. Low Pr means effective heat conduction with 
dominant thermal diffusivity. High Pr means effective heat 
convection with dominant momentum diffusivity. 

Ratio of inertial forces to viscous forces. Determines 
whether flow is laminar or trubulent. 


Table 9 

Definitions of some of the common dimensionless parameters employed in the analysis of phase change problems. 
Number Source Formula 
Biot, Bi [126] Bi =H! 
Nusselt, Nu [75,84] Nu =" 
Stefan, Ste [127] Ste = DT 
Dimensionless time or Fourier number, Fo [75] C= a 

p 

Rayleigh, Ra [84] Ra = LAN 
Prandtl, Pr [84] Pr=2 
Reynolds number, Re [122] Re=> 
Grashof number, Gr [128] Gr = gBATÉ 


v 


Approximates the ratio of buoyancy force to the viscous force. 


Table 10 


A list of some of the correlations derived by different references using dimensionless numbers. 


Reference(s) Correlation derived Interpretation of correlation 
84 V =a(t—te) Linear function of the melted volume of PCM with time. 
Vo 
É = 4,73F00-%%6$tel-538R 30-002 Melted volume fraction for the entire duration of test runs in 
à terms of Fo, Ste and Ra. 
Nu = 0.219 Ra0357 Pr0.019 (HH) 0067 Time averaged Nusselt number for the entire melting process. 
[74] ts = ro il + (0.25 + 0.17w°7)Ste} Total heat input per unit surface area Q and average heat flux 
at r=! over the duration of the melting process ty. 
q = KAT {1 + (0.1210 + 0.0424w)Ste(.7645-0.2022w) + 
1/4 
[78] keq(T) = 0.228 Ra(t)'/* [1 = a0) í Equivalent thermal conductivity and dimensionless 
melting rate [1 —A(T)]. 
A(T)[1 — InA(t)] = 1 — 4 fo Keq(t)dt 
[108] M(Ë) = LÉ fio (¥, Ddy, He(t) = ifs 2dy C(Ë) = iB æ dy with M=He —C Instantaneous melting rate, instantaneous heating rate and 
instantaneous cooling. 
= 2.39 x 10 “Fo®’® Ste® e? adial temperature gradient in the unfinned and finned tube 
97 aes 2.39 x 10 2F007655te0.528 Re0.195 Radial di in th finned and finned tub: 


OL = -3 F00.781 Ste0.616 Re0.503 
Tax = 2-23 x 10 ~Fo””*!Ste?-°!Re! 


PCM (MgCl-6H20) system. 


Ste were the important factors to consider in the phase change 
process. 

The difficulty in comparing data was illustrated by Agyenim 
[96] by comparing the dependence of Nusselt number on Reynolds 
number for four independent investigators: 


i. Hamada et al. [75] employed equation from Petukhov [129], 
Gnielinski [130]. 


(f jric/2)(Re — 1000) Pr 
14 12.7,/ f fic/2(Pr?? — 1) 
(3000 < Re < 10°, 0.5 < Pr < 2000)(1) 


Nu = 


ii. Zhang and Faghri [111] employed empirical correlation 
proposed by Edwards and Jensen [131]. 
102 
Nu = 0.023(1/d;) 


0.8 
=0350351/4) 0 (2) 


with l/d; calculated from available data as 0.50408. 
iii. Sexana et al. [132] calculated the Nusselt number using 
equation 3.34 for turbulent flow. 


Nu = 0.026Re°£ Pr! (3) 


iv. Esen et al. [80] and Rohsenow et al. [133] employed Dittus- 
Boelter equation for turbulent flow: 


Nu = 0.023 Re°®pr° for (Re > 2200) (4) 


Agyenim [96] illustrated that although all correlations showed the 
same general trend with an apparent slight differences in values 
recorded for the lower Reynolds numbers in the turbulent regime, 
major deviations increased with increased Nusselt and Reynolds 
numbers. Independent correlations presented demonstrated that in 
the turbulent region, the Nusselt number would be governed by the 
dimensionless parameters, Nu = Nu(Re, Pr) but data correlated only 
marginally well. 


5.1. Phase change problem formulation 


The analysis of heat transfer problems in phase change 
processes is complex due to the fact that the solid-liquid 
boundary moves depending on the speed at which the 
latent heat is absorbed or lost at the boundary. The position 
of the boundary is unknown and forms part of the solution Li et 
al. [134]. An increasing number of studies have been conducted 
and literatures on LHTES problems are numerous. Solutions to 
phase change problems include analytical, experimental and 
numerical using one-dimensional, two-dimensional or three- 
dimensional models to solve energy formulated equation. 
Energy equation is formulated in various ways with the phase 
change being accounted for in different representations. 
Examples of some of the formulation by various references 
are given in Table 11. Table 12 provides a list of references for 
the various solution methods employed in solving the phase 
change problems. 
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Table 11 
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Some examples of formulation for the phase change problem. 


Reference(s) Governing equation 
[74] T & 8 (rw 8) forr(t) <r < LT = Te, r < r(t) 
0 foraPCMslab insulated at one end 
where 1 +w = A/V with w = 4 1 foraPCM cylinder 
2 foraPCM sphere 

[s6 #12 (er) talm — pan 

pCpmr? À = 27rr;U(T — Ty) — mC, 2 

Phase change: H(T) = h(T) + p; f Ah; where h(T) = ifs pC pdT 

2 

[101] uC nas vite) Dane 

aH é AIN an. ar, 

C= OE (kor 7) ré (kp a) 
[108] Solid phase of PCM: # + uo % = œs ČT, liquid phase of PCM: + % =0 v? + gB(T — Ty) =0, © + À (uT) + (= a PT 
[121,135] Mass conservation equation: n =0 

Momentum conservation equations: o(% +u; a) i pgB(T — Tref) ee where oj; = — pô; + H (a H z) and p' + pgx2 

Conservation of energy: pc (@ +u; g) = Ea (x 2) 
[112] D E (kr D +18 (HD) 

Phase change formula: solid region (T<T,, —e): H=c,T, interface (Tn —Ee <T<Tn+e): H = CpT + ()T — Tm + €, liquid region (T>T,,+@): H=c,T+A 

Table 12 
Numerical, experimental and analytical solution and validation of phase change problems. 

Reference(s) Solution and validation Nature 
[66] 2D fixed-grid conjugate enthalpy based method. Finite difference and tridiagonal matrix algorithm solver. Numerical and experimental. 
[66] 1D enthalpy based method. Gauss-Seidel iteration procedure. Numerical and experimental. 
[80] 2D conjugate heat transfer enthalpy based method. Gauss-Seidel iteration procedure Numerical. 
[74] 1D enthalpy based method. Solved via Crank-Nicholson scheme. Numerical. 
[104] 1D enthalpy based method. Finite difference program-CYCLOP. Analytical and Numerical. 
[97] - Experimental. 
[93] 1D and 2D enthalpy method. Implicit finite difference method. Numerical and experimental. 
[102] 1-Parabolic and exponential approximations for temperature profile. Newton-Raphson technique. Analytical. 
[75] 2D enthalpy method. Solved using the control volume method by Pantakar (1980). Numerical and experimental. 
[105] 2D fixed-grid conjugate enthalpy based method. Control finite difference method. Numerical 
[81] 2D enthalpy based method. Gauss-Seidel iteration procedure and TRNSYS. Numerical 
[72] 2D Navier-Stokes and energy equations. Finite element model. Numerical and experimental. 
[68] 1D and 2D fixed-grid enthalpy based method. Control volume approach by Pantankar (1980). Numerical. 
[86] - Experimental. 
[91] 1D enthalpy based method. Experimental. 
[95] - Experimental. 
[76] Fixed-grid enthalpy-porosity model. Streamline upwind/Petrov Galerkin finite element method. Numerical. 
[106] 2D fixed-grid enthalpy method. Finite volume approach. Numerical 
[77] 2D enthalpy method. Solved using Ozisik. Theoretical 
[107] 1D enthalpy method. Analytical. 
[78] - Analytical and experimental. 
[108] 2D mass, momentum and energy equations. Finite difference method. Theoretical and experimental. 
[79] 1D and 2D enthalpy method. Fully implicit finite difference method. Numerical 
[101] 2D energy based equation. Finite element analysis. Numerical. 
[60] 2D fixed-grid enthalpy based method. Gauss-Seidel iteration method. Theoretical. 
[90] 3D enthalpy based method. Control volume method by Patankar. Experimental and numerical. 


A survey of phase change formulation reveals that the most 
common approach applied in the solution of phase change 
problems is the formulation using the enthalpy method. The 
advantage of employing the enthalpy procedure has been reported 
as the removal of the need to satisfy conditions at the phase change 
front which most authors indicated facilitate the implementation 
of the numerical algorithm. The survey of publications from LHTES 
systems with cylindrical containers also shows that nearly half 
(47%) of the studies undertaken are numerical and half of the 
references employed two-dimensional models (Fig. 4). 

Results from the numerical methods reported appear to show 
that they offer a good approach to solving the phase change problem 
although most of the available solutions to phase change problems 
apply to one- or two-dimensional systems due to the complexity of 
the equations involved in the phase change. The latent heat phase 


change problem has been assumed to be conduction controlled by 
some authors whereas natural convection in the liquid phase has 
been included in the analysis by other authors, adding to the 
complexity of the phase change problem. Few examples of phase 
change problem classification have been reported in Table 13. 
The numerical formulation has also concentrated on resolving 
specific phase change material problems and to study the 
characteristics of proposed or new materials. They include: 


e heat transfer during melting and solidification of PCM; evalua- 
tion of the temperature distribution inside PCM, interface 
transient location, and the effect of natural convection in the 
molten zone and, 

e the overall thermal performance, system heat capacity, charging 
and discharging rates and storage temperatures. 
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3-D na 
| Analytical 4% 14% 
Experimental 24% 
29% 
n.a — not 
available 


9 1-D 
50% 32% 


Numerical 
47% 


Fig. 4. Percentages of PCM solutions/analysis by methods and dimensions for 28 
LHTES systems with cylindrical containers. 


Table 14 


Table 13 
PCM heat transfer control classification by different references. 


Reference(s) PCM Control classification 
[68 Water-ice system Conduction 

[95 95% purity palmitic acid Conduction 

[86 95% purity lauric acid Conduction 

[76 n-Octadecane Convection 

[106] - Convection 

[7/71 n-Octadecane Conduction/convection 
[107] Paraffin wax, paraffin (p-116) Conduction/convection 

and stearic acid 
[78 n-Octadecane Conduction/convection 


Summary of method of analysis with key findings for conduction/convection control phase change problems. 


Reference(s) PCM/method Parameters used to examine results Key findings 
[86] Lauric acid. Pipe model with the PCM on the Reynolds and Stefan numbers. Graphical Charging time was less than solidification 
shell side and the heat transfer fluid flowing presentation of temperature variations in time due to: (i) melting rate being governed 
through the center of the tube. radial and axial directions for different by convective heat transfer and (ii) 
inlet HTF temperatures. conduction thermal resistance for the 
solidification process (Req = ry In(r/ry)/2ks). 
[78] n-Octadecane. PCM inside an isothermally Ra and Pr numbers. Photographic Melting proceeded faster than that of the pure 
heated horizontal cylinder. samples at dimensionless times, T. heat conduction and heat transfer at the upper 
Graphical presentation of equivalent region of the PCM was larger than the lower 
thermal conductivity equation derived to region due to the presence of natural 
estimate dimensionless melting. convection. However, decreasing Ra decreased 
convection and was negligible at Ra < 10°. 
[108] n-Octadecane. PCM enclosed in a rectangular Nu, thermal conductivity ratio A, solid A has practically no effect on the heating rate 


container and heated at a constant rate from 
one side. 


subcooling B and Ste, Heating rate He 
and Melting rate, M. 
Graphical presentation of Nu, A, B, M and He. 


He during the convection dominated regime 
and has a minor impact on He during the 
conduction dominated regime. Higher B values 
produce lower melting velocities, especially in 
the early stage of the process. Melting rate 
decreases as Ste increases. 


5.2. Classification of heat transfer characteristics 


Problems that arise due to heat transfer accompanied by a 
change of phase can be classified as: 


a. conduction controlled phase change, 
b. convection controlled phase change or 
c. conduction/convection controlled phase change. 


References to some of the control classifications for some of the 
various studies reviewed are given in Table 13. 

Studies based on the assumption that conduction is the major 
mechanism of heat transfer through the PCM assumed negligible 
convection heat transfer in the melt. In the case of convection 
controlled phase change, the common assumption was that the 
PCM was at the melting temperature. Models based on pure 
conduction and pure convection have not been able to predict 
accurately the melting rates nor track the complex motion of the 
solid-liquid interface at all times during the phase change. Less 
than 10% of articles surveyed treated convection and conduction at 
different stages and concluded that it depends on the phase 
transition stage to determine whether or not conduction or 
convection controls the heat transfer. Among these are Sari and 
Kaygusuz [86], Hirata and Nishida [78], Agyenim et al. [99] and 
Zhang and Bejan [108]. 

The authors accounted for both conduction and convection, 
making specifications at which times conduction/convection is 
(are) dominant during the phase change process. The details of the 
work by Zhang and Bejan [108] together with Sari and Kaygusuz 
[86], Agyenim et al. [99] and Hirata and Nishida [78] establish the 
following stages for the melting process: 


a. Sensible heat provision at the start-up period during which heat 
added increases the temperature of the solid PCM. Pure 
conduction heat transfer occurs at this stage. 

b. A second conduction regime, in which heat is purely transferred 
by conduction from the heated wall to the PCM and solid-liquid 
interface when melting had just begun. 

c. The transition from conduction to natural convection which 
starts as the thickness of the melt layer increases and the 
interface starts to incline. At the interface, there exists at 
equilibrium, a solid body and a pool of its own liquid. 

d. The convection regime, when most of the solid have been 
melted and the liquid core temperature distribution depends on 
height and depth and not on time (caused by buoyancy). 


Phase change problems, first treated as pure conduction controlled, 
has in recent times moved to a different level of complexity with 
added convection in the melt being accounted for (Table 14). 


6. Conclusion 


Studies of phase change systems have over the past three 
decades investigated design fundamentals, system and process 
optimization, transient behaviour, and field performance. The 
research and development has been broad based and productive, 
concentrating on both the resolution of specific phase change 
materials and problems and the study of the characteristics of new 
materials. The review reveals the different heat transfer solution 
methods employed by different researchers. They include theore- 
tical, experimental and numerical studies that have been 
conducted on the thermophysical properties of new phase change 
materials. 
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The following conclusions can be drawn from the models 
reported in the heat transfer and phase change problems reviewed: 


e Most of the research on phase change problems has been carried 
out within the temperature range 0-60 °C suitable for domestic 
heating/cooling. Common PCMs studied are paraffins, Glauber’s 
salts, wax, stearic acid and n-octadecane. Studies of medium 
temperature latent heat thermal energy systems to be applied for 
hot-side storage of absorption air conditioning are few. 

Most phase change materials have low thermal conductivities 
requiring heat transfer enhancement techniques to improve 
their rates of charging and discharging of energy. Majority of the 
studies undertaken on phase change materials have reported 
poor thermal conductivities requiring heat transfer enhance- 
ment techniques. As such current studies have concentrated on 
methods of heat transfer enhancement using mostly fins, 
insertion/dispersion of high thermal conductivity materials, 
multitubes and micro- or macro-encapsulation 

Numerical methods appear to offer a good approach to solving 
the phase change problem although most of the available 
solutions to phase change problems apply to one- or two- 
dimensional systems due to the complexity of the equations 
involved in the phase change. 

The latent heat phase change problem was first assumed to be 
conduction controlled. In recent times however, natural con- 
vection in the liquid phase has been included in the analysis of 
the phase change problem adding to the complexity of the phase 
change problem. 

e Acommon approach observed in the solution of phase change 
problems has been the use of the enthalpy formulation. The 
advantage of the enthalpy procedure has been reported as the 
removal of the need to satisfy conditions at the phase change 
front to facilitate the implementation of the numerical algo- 
rithm. 

Variations are found to exist in correlations derived to 
characterise the effect of melting and solidification of PCMs. 
The variations have been attributed to the different methods of 
analysis and the choice of characteristic lengths used and as such 
calls for caution in the selection of correlation for design 
purposes and quantitative comparisons. 

There are no national and international standards for testing 
PCMs. As such researchers have over the years applied different 
approaches in both test procedures and analysis of results. The 
development of a unified international standard for test and 
analysis will be helpful in the development of latent heat thermal 
energy storage systems. Also contradictions exist in the 
thermophysical properties of PCMs provided especially the 
latent heat of fusion, thermal conductivity and densities in solid 
and liquid states. This is again due to the absence of unified 
certification standards and procedures. 
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